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Abstract 

In this paper we give necessary and sufficient conditions on a nonnegative 
tensor to be diagonally equivalent to a tensor with prescribed slice sums. These 
conditions are variations of Bapat-Raghavan and Franklin-Lorenz conditions. 
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1 Introduction 

For a positive integer m let (m) be the set {1, ... , m}. For positive integers d, mi, . . . , rrid 
denote by ]]j m i x --- xm d the linear space d-mode tensors A = [<Hi,ia,...,id]i *j e ( m j)i3 € 
(d). Note that a 1-mode tensor is a vector, and a 2-mode tensor is a matrix. Assume 
that d > 2. For a fixed ik € (nik) the (d — l)-mode tensor [dii,...,i d ],^ € ( m j),j € 
{d)\{k} is called the (k,ik) slice of A. For d = 2 the slice and the (2,j) slice 
are the i — th row and the j — th column of a given matrix. Let 

s k,i k : = a iu-,ui ik e ( m k),ke (d) (1.1) 

i J &(m J ),je(d)\{k} 

be the (&,«'& )-slice sum. Denote 

Sk ■= (sk,i, ■ ■ ■ , Sk,m k ) T , k € (d) (1.2) 
Note that (fc,ifc)-slice sums satisfy the compatibility condi- 

mi m d 

J] a Mi = • • • = X] Sd ^- ( L3 ) 
il=l i d =l 

'This research started during author's participation in AIM workshop "Nonnegative Matrix 
Theory: Generalizations and Applications", December 1-5, 2008. 
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Two ci-mode tensors A = [aii,i 2 ,...,ij># = [bn,i a ,...,i d ] £ M m i x --- Xm d are called 
positive diagonally equivalent if there exist = (2^1, . . . ,Xk,m k 

)T G ^n k)k £ (d) 

such that Oi!,...,^ = bi 1 ,...,i d e Xl ' il+ "' +Xd ' id for all ij £ (nij) and j £ (d). Denote by 
j|mix...xm d ^ e cone of nonnegative,(entrywise), d-mode tensors. 

In this paper we assume that B = [&t 1) j 2 ,...,i <J ] £ j^™i x --- xm d i s a given nonnegative 
tensor with no zero slice (k,ii~). Let £ M™ fc ,£; € (d) are given k positive vectors 
satisfying the conditions fjl .3j> . Denote by R™ 1 x "' xmd (<?3, si, . . . , s^) the set of all 
nonnegative A = [aji,i 2 ,...,ij £ ]g)™i x --- m d h. av i n g the same zero pattern as B, i.e. 
&ix,...,i d = &ii,...,i d = for all indices ii, . . . , i^, and satisfying the condition 

(jl.ip . The aim of this paper is to give new necessary and sufficient conditions on B 
so that M™ lX '" m<i (S, si, . . . ,Sd) contains a tensor A, which is positively diagonally 
equivalent to B. For matrices, i.e. d = 2, this problem was solved by Menon [9] and 
Brualdi [I]. See also [ID]. For the special case of positive diagonal equivalence to 
doubly stochastic matrices see [5] and |13j . The result of Menon was extended for 
tensors independently by Bapat-Raghavan [2] and Franklin-Lorenz [6j. (See [I] and 
[12j for the special case where all the entries of B are positive. )0 In this paper we 
give a different necessary and sufficient conditions for the solution of this problem. 

Theorem 1.1 Let B = [bi 1: i 2 ,...,i d ] £ ]g>™i x --- xm d^ ^ > 2) ; be a given nonnegative 
tensor with no (k, ik)-zero slice. Let £ k = 1, . . . , d be given positive vectors 
satisfying &1.3\) . Then there exists a nonnegative tensor A £ jj™ lX --- xm d ? which is 
positive diagonally equivalent to B and having each (k,ik)-slice sum equal to Sk,i k , if 
and only the following conditions. The system of the inequalities and equalities for 
Xfc = (xfc,i, . . . ,x fcjmfe ) T eM. mk ,k= l,...,d, 

£i,u + x 2 ,i 2 + • • • + x d>id < if bi 1)i2) ... jid > 0, (1.4) 
sjxfc = for k = 1,... ,d, (1.5) 

imply one of the following equivalent conditions 

1. xi th + x 2 ,i 2 + • • • + x djid = if 6ti,i 2 ,...,id > 0. 

^- Z)6 4l ,i 2 ,...,j d >0 x Mi + X2 .*2 + • • • + ^d,i d = 0. 

Ln particular, there exists at most one tensor A £ jj™ lX --- xm d with (k,ik)-slice sum 
s k,i k for all k,ik, which is positive diagonally equivalent to B. 

The above yields the following corollary. 

Corollary 1.2 Let B £ j£™ lX --- xm ^ ^ > 2), fre a given nonnegative tensor with 
no (k, ik)-zero slice. Let s& € M™ fe , k = 1, . . . , d be given positive vectors. Then there 
exists a nonnegative tensor C £ flj™ 1 ' — m <^ which is positive diagonally equivalent to 
B and each (k,ik)-sum slice equal to Sfc,« fc > */ and only if there exists a nonnegative 
tensor A = [a-ii,i 2 ,...,i d ] £ ]g>™i x --- xm d^ having the same zero pattern as B, which 
satisfies 

For matrices, i.e. d = 2, the above corollary is due Menon [9j. For d = 3 this 
result is due to Thm 3] and for d > 3 [6j. Brualdi in [4] gave a nice and simple 

1 \ thank Ravi Bapat for pointing out these references in December 2009. 
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characterization for the set of nonnegative matrices, with prescribed zero pattern 
and with given positive row and column sums, to be not empty. It is an open 
problem to find an analog of Brualdi's results for d-mode tensors, where d > 3. 

Note that the conditions of Theorem 11.11 are stated as a linear programming 
problem. Hence the existence of a positive diagonally equivalent tensor A can be 
determined in polynomial time. If such A exists, we show that it can be found by 
computing the unique minimal point of certain strictly convex functions /. Hence, 
Newton method can be applied to find the unique minimal point of / and its value 
very fast. (See §3.) 

2 Proof of the main theorem 

Identify R mi x R m2 x . . . x R md with R n , where n = J2k=i m k- We view y G R n as 
a vector (xj~, . . . ,xJ) T , where x^ G R mk , k G (d). Let ||y|[ := \J y T y. Define 

/(y) = /((x7,...,xJ) T ) := K-,^ + - +Xd ^- (2-1) 

ije{m J ),je{d) 

Clearly, / is a convex function on R n . Denote by U(si, . . . , s^) C R n the subspace 
of vectors (pcj , . . . ,xJ) T satisfying the equalities (11.51) . 

Lemma 2.1 Let B = [bi lt i 2 .jJ G jg>™i x -" xm d^ ^ > 2), be a given nonnegative 
tensor with no (k, i k )-zero slice. Let s k G R™ k , k = 1, . . . , d be given positive vectors 
satisfying i fi.gj) . Then there exists a nonnegative tensor A G R™ 1 x --- xmd j which is 
positive diagonally equivalent to B and having each (k,i k )-slice sum equal to s k ^ k , 
if and only the restriction of f to the subspace U(si, . . . , s^), f/|U(si, . . . , s^), has 
a critical point. 

Proof. Assume first that /|U(si, . . . , s^) has a critical point. Use Lagrange 
multipliers, i.e. consider the function / — Ylt=i s fc x fc' *° deduce the existence of A = 
(Ai,...,A d ) T and (£,J , . . . , £,J) T G U(si, . . . , s d ), where £ k = . . . ,Ck,i k ) T , k G 
(d), satisfying the following conditions 

b iu ... tid e^ + - + ^i = \ k s k , ik , , i k G {m k ),ke (d). (2.2) 

je(d)\{k} 

Since > is a positive vector and the (k, i k )-slice of B is not a zero slice we deduce 
that X k > 0. Summing up the above equation on i k = 1, . . . , m k , and using the equal- 
ities CCl we deduce that X 1 = ... = X d >0. Then A = [b il ,... : i d e^ 1 ^- losXl)+ ^ + - + ^]. 

Vice versa suppose A = [bi 1: ...^ d e Xl ' il+ '" +Xd - l <i] has (k, i k )-slice sum equal to s k ^ k 
for all (k,i k ). Let l m = (1, . . . , 1) T G R m . Then there exists a unique ti G R such 
that s^T(xfc — t k l mk ) = for k G (d). Let £ fc := x^ — t k s k , k G (d). Then (|2.2p holds. 

□ 

Denote by V(si, . . . , s^) C U(si, . . . , s^) the subspace of all vectors (x^, . . . , xJ) T 
satisfying the conditionfTlof Theorem fTTTl Clearly, for each y G R n the function / has 
a constant value /(y) on the affine set y + V(si, . . . , s<j). Hence, if r] G U(si, . . . , s^) 
is a critical point of /|U(si, . . . , s d ) then any point in r} + V(si, . . . , s^) is also a crit- 
ical of /|U(si, . . . ,Sd). Denote by V(si, . . . ,Sd) C U(si, . . . , s^), the orthogonal 



3 



complement of V(si, . . . , s^) in U(si, . . . , s<j). Thus, to study the existence of the 
critical points of /|U(si, . . . , s^), it is enough to study the existence of the critical 
points of /|V(si, . . . , s c i) ± . Since the function e at is strictly convex for t G R for any 
a 7^ 0, more precisely (e at )" = a 2 e at > we deduce the following. 

Lemma 2.2 Let B = [^i 1 ,i 2 ,...,i ( i] € JR™ 1 ••• ' m <^ > 2), be a given nonnegative 
tensor with no (k, i\^)-zero slice. Let s*. G R™ fc , k = 1, . . . , d be given positive vectors 
satisfying hi. 3]) . Let U(si, . . . , s^), V(si, . . . , s^), V(si, . . . , s^) -1 be defined as above. 
Then /|V(si, . . . , s^)^ is strictly convex. More precisely, the Hessian matrix of 
/|V(si, . . . , Sfi) 1 - has positive eigenvalues at each point o/V(si, . . . , s^) -1 -. 

Theorem 2.3 Let B = \fii lt i 2r ., t i d \ £ jg>™i x --- xm d^ ^ > 2) } be a given nonnegative 
tensor with no (k, ik)-zero slice. Let Sf~ G R™ fc , k = 1, . . . , d be given positive vectors 
satisfying \1.3\) . Then the following conditions are equivalent. 

1. /|U(si, . . . , Sd) has a global minimum. 

2. /|U(si, . . . , Sd) has a critical point. 

3. lim/(y/) = oo for any sequence yi G V(si, . . . , s^)^ such that lim |jy/|| = oo. 
4- The only y = (x^, . . . , xJ) T G V(si, . . . , Sfc) that satisfies {l-4\ l is y = n . 

Proof. U El Trivial. 
[2]=>[3] Let (3 G U(si, . . . , s^) be a critical point of /|U(si, . . . , s^). Hence any point 
in (3 + V(si, . . . , Sd) is a critical point of /|U(si, . . . , s^). Hence /|U(si, . . . , s^) 
has a critical point £ G V(si, . . . , s^) ■ In particular, £ is a critical point of 
/|V(si,...,s d )- L . Let z G V(si,...,s d )- L ,||z|| = 1. For t G R define g z (t) := 
/(£ + tz). So g z is strictly convex on R and g' z (0) = 0. Let H(/)(y) be the Hessian 
matrix of /|V(si, . . . , Sd)' L at y G V(si, . . . , s^)^, i.e. the symmetric matrix of the 
second derivatives of /|V(si, . . . , Sd) ± at y G V(si, . . . , Sd) ± . Lemma 12.21 implies 
that the smallest eigenvalue a(y) of H(/)(y) is positive. Clearly, H(y) and hence 
a(y) are continuous on V(si, . . . ,Sd) ± . Hence miii|| y _^|| <1 a(y) = 2a > 0. There- 
fore, g'l(t) > a for t G [-1,1]. In particular g' z (t) > 2at and g z (t) > /(£) + at 2 
for any t G [0, 1]. So g z (l) > /(£) + a. Since g' z (t) increases on R it follows that 
g' z {t) > 2a for t > 1. Hence g z (t) > /(£) + a + 2a(t - 1) = /(£) + a(2t-l) for t > 1. 
Thus /(£ + u) > a(2||«|| - 1) for any u G V(si, . . . , s d ), ||u|| > 1. Hence E] holds. 
[3]=>[0 Since / = oo on dV(si, . . . , s^) 1 - it follows that /|V(si, . . . , Sd) 1 - achieves its 
minimum at £ G V(si, . . . , Sd) ± . Clearly, for any point y G U(si, . . . , Sd) there exits 
z G V(si,. . . ^rf)- 1 such that y G z + V(si, . . . , s d ). Recall that /(y) = /(z) > /(£). 
Hence /(£) is the minimum of /[U(si, . . . , s^). 

[3] [H Assume to the contrary that there exists / y = (x^, . . . , xJ) T G 
V(si, . . . , Sd)- 1 - which satisfies ()1.4|) . Hence, there exists ij G (rrij),j G (d) such 
that feij,...,^ > and + ... + x d ,i d < 0. Thus, there exist ai,...,a p < 
and ... ,Pi > such that /(ty) = 7 + X)f=i Pi etai ■ (Each a q is equal to some 
+ • • • + Xd,i d < 0, where fei lv ..,i d > 0, and each /3 g is a sum of corresponding 
bi lt ...,i d > 0.) Hence, limt_ >tX) /(iy) = 7, which contradicts [3] 
[1 4'|31 Let y G (xj , . . . , xJ) T G V( Sl , . . . , s^, \\y\\ = 1. Then 

7i(y) := max a?i ^ + . . . + x di i d > 0. 
6i 1 ,...,i d >0 
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The continuity of h(y) on the unit sphere in V(si, . . . , s c i)' L implies that 

min My) = a > 0. 

yeV(s 1 ,...,s d )M|yj|=l 

Let f3 = niiiu, ; . {) b ; . > 0. Hence, for any y £ (xj , . . . ,xJ) T € V(si,. . ^Sd) 1 

||y|| = 1 and t > we have that f(ty) > f3e at . This inequality yields [3j □ 

Proof of Theorem II. 1L Assume first that there exists a nonnegative tensor 
A € ]^™ lX --- xmd ) which is positive diagonally equivalent to B and having each (k, i^)- 
slice sum equal to Sk,i k - Lemma I2TT1 yields that /|U(si, . . . , s^) has a critical point, 
i.e. the condition of Theorem 12.31 holds . Since /|V(si, . . . ,s ( j,) ± is strictly convex, 
it has a unique critical point £ € V(si, . . . , s ( i)~ L - Hence all critical points of a 
convex /[U(si, . . . , s^) must be of the form £ + V(si, . . . , s^). The proof of Lemma 
12.11 yields that A is unique. 

Theorem 12.31 implies the condition [7} Hence the conditions (jl.4p and (|1.5p yield 
the conditions U\ and of Theorem 11.11 

Assume that the conditions (|1.4p and (jl.5p hold. Clearly the conditions [7] and 
[Hof Theorem 11.11 are equivalent. Suppose now that the conditions (II. 4p and (11.5P 
imply the condition [I] of Theorem ll.il Hence the condition [^] of Theorem 12.31 holds. 
Use the the condition [D of Theorem 12.31 and Lemma 12.11 to deduce the existence of 
a nonnegative tensor A £ R^ 1 "" m<i , which is positive diagonally equivalent to B 
and having each (A;,ifc)-slice sum equal to s^^ k . 

□ 

Proof of Corollary 11.21 We prove the nontrivial part of the corollary. Suppose 
that there exists a nonnegative tensor A = [a>i lt i 2 ,...,i d ] € flj" 11 *— h av i n g th e sam e 
zero pattern as B, which satisfies (II. ip . Clearly, A is positively diagonally equivalent 
to A and has each (k, ifc)-sum slice equal to Apply Theorem 1 1.1 1 to A to deduce 
that the set of inequalities +£2,i 2 + • • - + Xd,i d < if ai 1: i 2 ,...,i d > 0, together with 
the equalities (jl.5p yields the condition ^ . . >0 ^i,ii + £2,i 2 + . . . + Xdi d = 0. 
Since a^...^ > <^=^ hi,...,^ > we deduce that the conditions (jl.4j) and ([1.5p 
of Theorem 11.11 yield the condition of Theorem 11.11 Hence there exists a nonneg- 
ative tensor C € ]^™ lX --- xm <i ) which is positive diagonally equivalent to B and has 
(fc, ifc)-sum slices equal to Sfc,j fe . □ 



3 Remarks 

Theorem 11.11 the main result of this paper, is stated stated in terms of linear pro- 
gramming. Hence by the results of [H [7] one can verify these conditions in poly- 
nomial time. The proof of Theorem 12.31 combined Lemma 12.11 shows that to find 
A, which is diagonally equivalent to B, we need to find the minimum of the strict 
convex function /|V(si, . . . ,s ( i)' L - There are many numerical methods to to find 
the unique minimum, e.g. |3j. Since the Hessian at the critical point of our strict 
convex function has positive eigenvalues, one should use the Newton algorithm, or 
its variant as Armijo rule [11J, to obtain the quadratic convergence. 

In the special case of diagonal equivalence to doubly stochastic matrices, one can 
performs the Sinkhorn scaling algorithm |13|. IS], which converges linearly. It seems 
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to the author, that even in the case of matrices, a variant of the Newton algorithm 
should outperform the Sinkhorn scaling algorithm. The numerical aspects of the 
comparison between the two algorithms will be done somewhere else. 
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